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Abstract 

I We consider black holes localized on the brane in the Randall-Sundrum infinite braneworld model. 

^S) ■ These configurations are static and charged with respect to a spherically symmetric, electric Maxwell 

field living on the brane. We start by attempting to construct vacuum black holes, in which case our 
conclusions are in agreement with those of Yoshino in JHEP 0901:068, 2009 ( arXiv:0812.0465 ). Ahhough 
approximate solutions appear to exist for sufficiently small brane tension, these are likely only numerical 
artifacts. The qualitative features of the configurations in the presence of a brane U(l) electric field are 
similar to those in the vacuum case. In particular, we find a systematic unnatural behaviour of the metric 
functions in the asymptotic region in the vicinity of the AdS horizon. Our results are most naturally 
interpreted as evidence for the nonexistence of static, nonextremal charged black holes on the brane. In 
contrast, extremal black holes are more likely to exist on the brane. We determine their near- horizon 
^ ' form by employing both analytical and numerical methods. For any bulk dimension d > 4, we find good 

I agreement between the properties of large extremal black holes and the predictions of general relativity, 

J I ■ with calculable subleading corrections. 

! 1 Introduction 

O 

The Randall-Sundrum (RS) infinite brane"W"orld scenario [1] has been proposed as a mechanism to explain 
the hierarchy bet"ween the TeV scale and the Planck scale, and to realize four-dimensional gravity effectively 
on the 3-brane. In this model, the observable universe is a 3-brane (domain "wall) to "which standard model 
' fields are confined, while gravity can access the extra spatial dimensions. The bulk metric is a locally Anti-de 

^ , Sitter (AdS) spacetime satisfying the Einstein equations "with negative cosmological constant. At lo'w enough 

energy, perturbative Ne"wtonian gravity is recovered on the brane at distances large compared to the AdS 
length scale £. 

The existence of black hole solutions in the RS model is an interesting open problem (see e.g. Kanti [5] 
for a recent revie"w on this issue). A priori, it is not clear why such black hole solutions should not exist, 
and the existence of bulk black hole solutions would imply that there are also black hole solutions on the 
brane. Indeed, such exact black hole solutions were found by Emparan et al. [3] for the dimension d = 4 
of the bulk spacetime, where the static black hole is localized on a 2-brane. The construction of such black 
hole solutions in [3] relied on the existence of a special class of solutions of the d — A Einstein gravity - the 
so-called C-metric |3], and the brane is a suitable slice in this spacetime. 

However, in the absence of a d > 4 generalization of the C-metric, the existence of higher dimensional 
brane world black holes is controversial. Although Dadhich et al. [S] proposed an analytic solution for brane 
black holes in the d = 5 RS scenario, this result was obtained subject to approximations whose validity is 
not clear. In that approach, the braneworld Einstein equations [6] were solved assuming a particular form 
for the bulk Weyl tensor on the wall. Interestingly, this resulted in a Reissner-Nordstrom (RN) geometry, 
although there is no brane electric charge. However, in this case only the induced metric on the brane was 
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found, without solving the bulk equations of motion. Thus the solution proposed in [S] for a brane black 
hole in the RS scenario is not really satisfactory. 

Hence, for bulk dimensions c? > 4, one has to rely on numerical methods to address the existence of black 
holes in the RS model. So far, the only results supporting the existence of such solutions on a nonperturbative 
level are those by Kudoh et al. [71 [S]. They were obtained within an in principle rigorous approach, i.e. 
by solving the full set of equations in the bulk with suitable boundary conditions on the brane. Reporting 
numerical results for spacetime dimensions d — 5 and 6, Kudoh et al. [71 [B] presented localized black holes 
whose horizon radius was much smaller than the AdS length scale. In contrast, for larger black holes their 
numerical accuracy decreased and convergence was lost within their numerical scheme. 

Moreover, the recent numerical work of Yoshino 9_ (while finding agreement with earlier perturbative [10] 
and numerical [3 [5] results) also failed to find such large black hole configurations. Performing a careful error 
analysis in terms of systematic and nonsystematic errors, Yoshino 9 concluded that the positive results in 
[3 [5] were nothing but numerical artifacts, mainly due to an inappropriate treatment of the outer boundary 
of the integration region. 

In fact, Emparan et al. [11) and also Tanaka '12| had earlier put forward a number of theoretical arguments 
against the existence of d > 4 static black holes on the brane in the RS model, which were mainly based on 
a version of the AdS/CFT correspondence. According to the conjecture by Emparan et al. [TT], such bulk 
black holes would necessarily be time-dependent, since their duals would describe quantum corrected black 
holes in a d — 1 dimensional braneworld. Counter-arguments were, however, given by Fitzpatrick et al. ^13) . 
Further results on black holes localized on the brane can be found in [T4]-|17j. 

Motivated by these confiicting arguments and results, we have performed an independent investigation 
of the issue of the numerical construction of d > 4 braneworld black holes. Similar to the previous work by 
Kudoh et al. [71IH] and Yoshino [S], we have solved the bulk Einstein equations with Israel junction conditions 
on the brane, employing different numerical methods from those in [71 |S1 Hj- In particular, we have used 
a compactified radial coordinate which in principle avoids the problems associated with the position of the 
outer boundary of integration for the radial coordinate. 

The results we have found for bulk dimensions d — 5 and 6 support the claim of Yoshino '9| that 'a 
solution sequence of a static black hole on an asymptotically flat brane that is reduced to the Schwarzschild 
black hole in the zero tension limit is unlikely to exist'. In particular, we have noticed a systematic unnatural 
behaviour of the metric functions for large values of the radial coordinate, i.e. close to the AdS horizon. 
This behaviour seems to be at the origin of the loss of numerical convergence for large black holes. 

Obviously, it is interesting to examine how generic this behaviour is, and whether it can be circumvented 
in more general cases. Perhaps the simplest more general case to be studied corresponds to black hole 
solutions, which are electrically charged with respect to a Maxwell field living on the brane. This type of 
solutions was considered by Chamblin et al. [TS], following a different approach, however. There the "initial" 
data on the brane was prescribed, and then it was evolved in the spacelike direction transverse to the brane, 
by solving the bulk equations numerically. The results in [18j show the occurrence of pathological features 
in the bulk for any initial data. However, Chamblin et al. [18] postulated a special restricted form of the 
braneworld metric, with a single essential function; moreover, the numerical integration employed a relatively 
small cutoff radius. 

Our present approach is rather different from 18 , since we attempt to directly solve the bulk vacuum 
Einstein equations with suitable boundary conditions. The bulk theory is the same as for the uncharged 
case, the electric charge entering the problem via the Israel junction conditions on the brane. Restricting to 
static, nonextremal configurations in an AdSs bulk, our results show that all pathologies present for vacuum 
black holes occur also in this case. Therefore, we conclude that the existence of static, charged, nonextremal 
solutions is unlikely, as well. 

Turning next to static, extremal, electrically charged black holes in the RS brane world scenario, how- 
ever, we anticipate that such solutions are more likely to exist, since a number of arguments put forward 
against the existence of static black holes do not apply when the Hawking temperature vanishes. Here we 
investigate the near-horizon structure of extremal black hole solutions with electric charge on the brane, 
without attempting to construct the full configurations. This restriction leads to a system of coupled non- 
linear ordinary differential equations, which are solved numerically within a nonperturbative approach for 
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several dimensions d > 5 of the bulk. 

For a five-dimensional AdS bulk, this problem has been considered by Kaus and Reall |19| . In this work 
we generalize their results for any d > 5 dimensions, and show that for large black holes, there is a good 
agreement for the induced metric on the brane and for the entropy with the predictions of general relativity 
(GR), with calculable subleading corrections. 

The paper is organized as follows: in the next Section we present our results for nonextremal black holes 
in the RS brane world model, employing a nonperturbative approach, by directly solving a set of three 
nonlinear partial differential equations with suitable boundary conditions. In Section 3 we consider extremal 
black holes that are charged with respect to a purely electric Maxwell field on the brane and determine 
their near-horizon form. We give our conclusions and remarks in the final Section. The Appendix contains 
a discussion of some technical aspects involved in our numerical investigation of nonextremal brane world 
black holes. 



2 Nonextremal configurations 
2.1 The problem 

We consider the RS braneworld model, with a d-dimensional bulk spacetime and a single {d— l)-dimensional 
brane with positive tension in it. Also, we impose Z2-symmetry about the brane, which is assumed to be 
asymptotically flat. The bulk matter is merely a negative cosmological constant, and the brane tension and 
the matter localized on the brane are treated in a distributional sense. 
The action of this model is 



WTrGdJM Jbranc \8nGd IGnGd-i J 



(2.1) 



where Ai is the bulk spacetime, A = — (d — 2)/{d — 1)/{2P) is the bulk cosmological constant and K^^ is 
the projection of the extrinsic curvature of the brane hyper-surface with induced metric Also Gd is 

Newton's constant in d-spacetime dimensions; a and F — dA are the brane tension and the field strength of 
the Maxwell field on the brane, respectively. 

From the above action, we obtain the d-dimensional Einstein equation in the bulk 

-Ry - l^mj + ^9^3 = 0. (2.2) 
For the RS infinite braneworld scenario, the Israel junction conditions on the brane are given by 

i^M^-i^V = 47rGd(-CTV + T-7^ V ) > (2-3) 

\ 47rGd-i / 

where t^^ is the energy-momentum tensor of the matter fields on the brane. For a U(l) field, the expression 
of t^i, is 

+ _p p " _ i p „ pa/3;, C2 A] 

We shall set the brane tension to the RS value a — while Gd-i — ^^Gd- This simplifies eq. (|2.3|) 

to 

d-2. 2i 
d~3 



Kh^i, — + ^tpi^n, (2-5) 



which is the form used in what follows. The brane U(l) field is a solution of the Maxwell equations 

Vi^F^" = 0, (2.6) 

for a metric background given by /i^i,. 
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2.1.1 The ansatz and the equations 

The metric ansatz employed here essentiaUy corresponds to the one used in the previous studies [71 [U |2] ■ 
The black hole metric is spherically symmetric on the brane and axisymmetric in the bulk spacetime, with 
line element 

ds^ = -^1-. (^e^Bir^^){^ + r^dx^) + e2C(-,x)^2 ^^^2 ^^^2_^ _ e^Air.x)F{r)dt^^ , (2.7) 

which is parametrized in terms of two background functions 

rQ\<^^^ r 



F(r) = 1 - (^-^i j , zir,x)^l + -^cosx, (2.8) 

and three unknown metric functions A{r,x), B{r,x) and C(r, x)- In the above relations, rp is a positive 
constant and dfi^ g is the metric on the {d — 3)-sphere. The background functions F{r) and z{r,x) have 
been introduced such that two important limits of the general solution are already contained within the 
ansatz (|2.7p . For I co one finds the well-known Schwarzschild-Tangherlini black hole, expressed in the 
usual Schwarzschild coordinatcfQ- Another limit of interest is tq = 0, in which case one recovers the original 
RS model, i.e. a part of AdS^ spacetime expressed in Poincare coordinates (and ^ = _B = C = Oin both 
limits) . 

The event horizon is supposed to reside at a surface of constant radial coordinate r — and characterized 
by the condition F{ro) — 0, while the brane is located at x = 7r/2. Then the coordinate range considered is 
To < r < oo and < x < 7r/2. Therefore the coordinates in eq. (|2.7p have a rectangular boundary and thus 
are suitable for the numerical methods employed. 

The induced metric on the brane has line element^ 



da^ = grr{r)dr^ + 5no(r)df7^_3 + gtt{r)dt^ , (2.9) 



with 



g2B(r,77/2) 



F{r) 

i.e. it describes a static, spherically symmetric black hole spacetime in d — 1 dimensions. 

The equations satisfied by the functions A, B and C are found by using a suitable combination of the 
Einstein equations, Gl + A — 0, + + 2A = and + A = (where (f> denotes an angle of the d — 3 
dimensional sphere) and reacd 

All 1 V ,d-2 iF' cosx, (d-3)F' cosx, , cosXso 

cosY,,^, cosy, (d— 3)cotx,; rsinv, 1 • rsinv.o 
Hd - 3)iA' - ^)(G' - ^) + ,4 \ a +^) + ^iA+ (2.11) 

(d--3) • rsinX wA rsinX s _ (d - l)e^^ 1 (r + £cosx) cos^ X ^ ^ 
r^F ^ Iz Iz ' eFz^ rF e^z^ i^z^ 



'^Refs. [T] Is] and [9] preferred to describe the Schwarzschild black hole in isotropic coordinates, which was also our initial 
choice. However, we realized that a Schwarzschild coordinate system together with the coordinate transformation (12.161 1 
improved the quality of the numerical calculations for d = 5. 

^Note the presence of three distinct metric functions in eq. 1 12.91 1. In contrast, in GR only two metric functions are present, 
with the usual choice grr = 1/Af(r), gnn = and gtt = N{r)cr'^(r). 

^One can see that the case d = 4 is special, since a number of terms vanish in this case. However, the structure of the 
equations is the same for any d > A. 
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where a prime denotes the derivative with respect to the radial variable r and a dot denotes the derivative 
with respect to the angular variable x- The remaining equations =0, G^ — = yield two constraints. 

Following gl], we note that setting Gj+A = 0, G;; + Gx + 2A = 0, G^+A = in V^G^'' = and V^G^^ = 0, 
we obtain Cauchy-Riemann relations for G^ and GJ! — G^. Thus the weighted constraints satisfy Laplace 
equations, and the constraints are fulfilled, when one of them is satisfied on the boundary and the other at 
a single point [H]. 

For completeness, we here present also the expressions for the Hawking temperature Th and the event 
horizon area of a bulk black hole, 

Th = e^(0'X)-^(«^^)^^^, A^'^ ^ Vd-.r^' T^' eS(o^x)+(rf-3)C(o,x) ^j^Llj^dx, (2.14) 

47r^o Jq (i + nicosx) 

where Vd-3 is the area of the unit 5'''"'^ sphere. (The Einstein equation G^ = implies that e"^~^ is indeed 
constant on the horizon.) The associated black hole on the brane would have the same Hawking temperature 
as the bulk solution, while its event horizon area would be 

A'-t'^ - T/d_3r^'e('^-3)c(o,V2). (2.15) 
In practice, we have found it convenient to introduce the radial coordinatcQ 



P=^Jr'-rl (2.16) 

such that the horizon would reside at p = 0, which gives a simpler set of boundary conditions for the black 
hole horizon. The transformation then results in a new form of the general ansatz (12.71) with 



^ (''''^'''^^:^ + + + e2^(^^'^)(p^ + r^) sin^ xd^U " e'^^^^^^ F,{p)dt''^ , (2.17) 

where Fi{p) = F{r{p)){p^ + rl)/p^ and F2{p) = F{r{p)). (Note that Fi{p) (d - 3)/2 + 0{p^), F2{p) ^ 
{d — 3)p^/(ro) + 0(p'') as p — >■ 0.) The temperature and the horizon area stiU given by p.l4p . (|2.15p . 



*This change of the radial coordinate proved useful before in the numerical study of non-uniform black string solutions I22| . 
which is an axisymmetric problem with some similarities to the problem studied here. 
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2.1.2 U(l) field on the brane and boundary conditions 

Perhaps the simplest example of nonvacumxi solutions in the RS infinite braneworld scenario is provided by 
a Maxwell field confined to the brane. In this work we shall restrict to a static, spherically symmetric, purely 
electric field with U(l) potential 

A = V{p)dt. (2.18) 

Thus the field strength tensor is F = ^dp A dt. The Maxwell equations (I2.6P imply the existence of the 
first integral 

Vip) = Q / dp ^__e^(PW2)+B(P.V2)-(d-3)c(p,V2)^ (2.19) 

where Q is an integration constant that fixes the electric charge on the brane. 

The numerical solution of the equations is pursued subject to the following set of boundary conditions 

on the black hole horizon, 

^1 =B\ =C\ =0, (2.21) 

at infinity, and 

= dxB\^=, = dxC\^^, = 0- (2-22) 

on the symmetry axis. Regularity at x = further requires that i?|x=o = C\x={)- This condition can be 



implemented by working with the new function C = C — B. The Israel junction conditions ()2.5p together 
with the expression (|2.19p for the U(l) potential lead to the following boundary conditions on the brane 



V^P'+rp ( B(p.tt/2) _ 1^ I 2d-7 B( p,7r /2)-2(d-3)C ( p.tt /2) I 

I ^) ^ (d-2){d-3)^ (p2+r2)(2d-7). 



X''^\x=7r/2 ~ '^X-^\x=7r/2 ~ I ■ ' ' (d-2){d-3)^ (p2+r2)(2d-7)/2 

Q _ y/p^ + rg ,B(p,^/2) 3 B(p,^/2)-2(d~3)C(p,7,/2) ^ ^2 /r, r,o^ 

Wx=.l2- I ^) (d-2)(d-3) (p2+^2)(2<i-7)/2^ ' i^-^^^ 

where Q = corresponds to the vacuum limit. 
2.2 Numerical results 

The results exhibited in this Section concern the physically most interesting case c? = 5 (i.e. a 3-brane). 
However, we have observed a similar behaviour when considering vacuum black holes for a 4-brane instead. 

The numerical calculations have been based on the Newton-Raphson method and been performed with 
help of the program FIDISOL/CADSOL [23], which also provides an error estimate for each unknown 
function. Different from previous work on this problem, in our approach the boundary conditions p.2ip are 
really imposed at p = c». This is achieved by employing a compactified radial coordinate x — p/il -\- p) 
which maps p = oo to the finite value x = 1. Further details on the numerical method are given in the 
Appendix. 

The first relevant input parameter for our problem is the dimensionless ratio 

L=j- (2.24) 

Without any loss of generality, by using a simple rescaling of p, one can set rp to take a fixed value, and then 
vary the AdS length scale £. The results reported in this Section correspond to the choice tq = 1, although 
similar results have been found for other values of rg. In the charged case, there is a second input parameter 
of the problem, which is the ratio 

Q 



'^0 



(2.25) 
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2.2.1 The vacuum case {Q — 0) 

Similar to Kudoh et al. [71 [S] and Yoshino |5] , we have first applied the numerical procedure to the above 
stated set of equations and boundary conditions, choosing a large value of the AdS length scale t (typically 
L ~ 10~* — 10~^) and employing the initial guesfrl A = _B = C = 0. In this case, the solver has converged 
and provided numerical output with good accuracjo- The functions A, B and C have a nontrivial shape, but 
their magnitude is small, almost zero. Next, we have started a new iteration and employed this data as initial 
guess for a slightly larger value of L. Again, the solver has converged and provided apparently reasonable 
numerical output, which in turn has been used as input for the next iteration for another slightly larger 
value of L, etc. For all configurations with sufficiently small values of L, we have found good agreement of 
the extracted physical properties with those reported by Kudoh et al. [3 18] , obtained within their numerical 
scheme. 

In principle, working in small steps, one expects to find solutions for arbitrarily large values of L in 
this way. However, for any grid choice employed, we have noticed that the numerical accuracy appeared 
to be progressively deteriorating with increasing L. Around L ~ 0.3 finally, the numerical errors have 
turned unacceptably large for the configurationf|j, while for still larger values of L the numerical errors have 
accumulated further, until finally convergence of the numerical scheme has been lost. 

In our scheme, the numerical problems appear to originate mainly in the asymptotic regior0. As L 
increases, the metric functions start to develop an increasingly unnatural shape for large values of the radial 
coordinate (typically for p > \Qt} and all values of x (although in some functions this feature is more 
pronounced for x ^ 0). This unnatural behaviour manifests itself in the occurrence of "oscillations" of the 
metric functions in the far fielc^. The amplitude of these "oscillations" increases with increasing L. At the 
same time these "oscillations" start at increasingly smaller values of the radial coordinate. We illustrate this 
unnatural behaviour of the metric functions close to the AdS horizon for a typical d — b calculation (with 
large numerical errors) in Figure 1 (left) and also in Figure 2 (the case with q — Q). 

This unnatural behaviour appears for all grid choices and all metric parametrizations considered, includ- 
ing the one employed in [TllHJI^. Also, quite interesting, when working with a finite grid with cut-off radius 
Pmax — lOr-Q (the typical case considered in [7]), one usually fails to detect this oscillatory behaviour, since 
it appears only for larger values of the radial coordinate. (A distinct dependence of the numerical results 
on the value of Pmax was already noticed by Yoshino , representing an essential part of the nonsystematic 
error present in his calculations.) 

Moreover, since the numerical error increases gradually, it is impossible to identify a critical value Lc, 
such that one could claim non-existence of the solutions for L > Lc- Thus the "oscillations" are likely to 
exist for any value of L. We surmise that for very small values of L such a behaviour is located at values of 
X very close to unity {i.e. at very large values of p), and that the amplitude of the "oscillations" is at the 
same time very small (and thus difficult to detect within the numerical approach employed). 

2.2.2 U{1) field on the brane 

It is interesting to examine whether the above results hold also for more general cases. In particular, one 
would like to know whether the unnatural asymptotic behaviour, that we have noticed in the vacuum case, 
survives in the presence of matter fields on the brane or whether it can be circumvented due to their presence. 

The simplest case one may think of corresponds to a spherically symmetric static Maxwell field living 
on the brane. As one can see from eqs. ()2.23p . although the bulk equations are the same as in the vacuum 
case, the presence of the U(l) field leads to a different set of boundary conditions on the brane and thus to 
a different geometry in the bulk. Naively, one may then expect the metric on the brane to correspond to a 
RN black hole with corrections. 

^We note that for Lj^O, A = B = C = Ois not a solution of the field equations unless ro = 0. 
®The computed relative error of the "solution" (truncation error) was on the order of 0.001 in this case. 
'^This has also been manifest in the constraint equations and Gp — G^, implying large errors in the evaluation of a 
Hawking temperature. 

*In previous numerical work numerical problems were largely associated with the symmetry axis x = 0. 
similar behaviour was noticed by Yoshino |9] (see his Figure 3). 
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Figure 1. The metric functions A, B, and C are shown as functions of the compactified radial coordinate 
X — p/{l + p) and the angular variable x for ^ vacuum configuration with L = 1 (left) and a charged one with 
L = 0.24, q = 1.73 (right). 

Our initial hope has been that the presence of an electric charge on the brane would make the numerical 
scheme more stable. However, in contrast to our expectations, our results for typical charged configurations 



8 



have turned out to be qualitatively similar to the vacuum casJ^. For a given value of L, we have employed 
the corresponding vacuum data as the initial guess for the numerical integration of the equations in the 
presence of a small electric charge Q. Subsequently, we have slowly increased the value of Q, and thus the 
second dimensionless parameter q. 

For small values of i (typically L < 0.1), the solver has converged and has provided apparently reasonable 
numerical output, as in the vacuum case. In particular, the shape of the functions appears to be rather 
insensitive to the value of q, while their magnitude is increasing with q. When extracting the Hawking 
temperature and the horizon area from the numerical data for these small values of L, we have observed 
that for fixed L the Hawking temperature decreases with increasing q as expected, while the horizon area 
increases for the "supposed" bulk black hole as well as for the associated "supposed" black hole on the brane. 

However, for any value of q employed, we have seen the numerical accuracy of the calculations to strongly 
deteriorate with increasing L as in the vacuum case, until at some stage the numerical solver has stopped to 
converge. Moreover, the behaviour of the metric functions for large values of p is also qualitatively similar 
to what we have found for q — 0, while the magnitude of the "oscillations" in the asymptotic region even 
increases with increasing q, as seen in Figure 2. This unnatural behaviour appears to be generic for any 
choice of the nonequidistant grid employed in the integration, and we conclude that all pathologies observed 
for vacuum configurations are also present in the charged case. Clearly, this puts the existence of such 
charged black hole solutions into strong doubt. 

We would like to emphasize that all of the configurations obtained numerically yield a nonvanishing 
Hawking temperature, and thus would correspond to non-extremal black holes (if they were indeed solutions) . 
Evaluated according to eq. (|2.14p . the expression for the Hawking temperature does not yield a strictly 
constant value, but typically varies slightly (because of the numerical error), as long as L and q are sufhciently 
small. For the charged configurations obtained here, however, the numerical error makes |A(0,x) ~ B{0,x)\ 
(and thus the temperature) deviate appreciably from a constant value. In fact, the variations become the 
larger the larger the charge q, whereas in principle the equations should guarantee a constant value of the 
temperature. It is thus not surprising to see the numerical integration fail to converge for some critical value 
of the charge. 

On the other hand, with increasing q, the "supposed" Hawking temperature of the configurations de- 
creases, and configurations with appropriately large values of q would be expected to correspond to near- 
extremal solutions. Indeed, this reasoning suggests that, similar to GR, extremal solutions might exist, which 
possess a maximal value of the electric charge for a given event horizon radius, and whose temperature would 
be zero. However, the study of black holes which are extremal or close to extremality is a difficult numeri- 
cal problem already for situations, where the existence of solutions is not controversial. Consequently, the 
existence of extremal braneworld black holes in the RS model is an open problem, which we cannot attack 
within our current numerical scheme. However, we will address the problem from a different direction in the 
next Section. 

2.2.3 The interpretation of the results 

We think that the interpretation of our results can only be analogous to the interpretation proposed by 
Yoshino j^. In particular, we would like to consider the following possibilities: 

First, one might suspect that our approach and/or the numerical methods are wrong. However, similar 
methodq^ were used in the past to solve a variety of problems, see e.g. [22l [24l [25]. Moreover, we have 
extensively tested these numerical routines to recover numerous exact solutions in GR and field theory. At 
the same time, some of the new solutions derived by using the code FIDISOL/CADSOL were rederived 
subsequently by other groups with different numerical methods. Therefore we think that this hypothesis can 
be safely excluded. 

The second possibility would be that static solutions exist, but that they are very hard to find, forming 
isolated "islands" in the parameter space. Moreover, these "islands" would be disconnected from the L = 

'^''We have considered charged configurations only in d = 5 dimensions. 

^^In particular, the treatment of the behaviour of the metric functions on the axis X = ^^'^ ^-t infinity has been similar to 
that used in this work. 
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Figure 2. The metric functions A, B and C are shown at x = (on the symmetry axis) and sd, x ~ '""/S 
(on the brane) for configurations for L = 10 and three values of the charge parameter q. 



Schwarzschild black hole used as initial guess in the numerical attempts to construct such solutionq:!^. In 



^^Such examples are known for numerical solutions in field theory, the case of vortons (which are d = 4 flat space toroidal 
solitons somewhat analogous to d = 5 black rings) being perhaps the most notorious. There several disconnected branches of 
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this setup, the numerical results found would be numerical artifacts, since the true solutions would require a 
starting profile different from the Schwarzschild black hole. However, we consider this possibility as counter- 
intuitive and unlikely. 

In our opinion, the most likely scenario is that there are no static solutions of the bulk Einstein equations 
(I^TT|) - (PTTO)) with the boundary conditions ([1201) -(12221) and ([2231) • Moreover, in agreement with the con- 
jecture put forward by Emparan et al. [11 and Tanaka [12^, all nonextremal black holes on the brane would 
be dynamical objects. If this were indeed the case, then the numerical results could be interpreted as fol- 
lows: what we have found should correspond to static approximations of some yet unknown time-dependent 
configurations. The full solutions for black holes on the brane would be dynamical, quantum corrected, 
evaporating black holes; thus also the bulk solutions would depend on time. The systematic occurrence of 
the observed unnatural behaviour of the metric functions close to the AdS horizon would be the result of 
an inappropriate static metric ansatz. For black hole with sizes small compared to the AdS length scale, 
the pathological behaviour at infinity would not be seen due to the lack of numerical accuracy. In other 
words, for small values of |A|, the hyperbolic character of the equations would not become manifest, and the 
solver would fail to see the quantum corrections. However, as the value oi L = r^jl increases, the neglected 
dynamical terms would destabilize the numeric£p^. This interpretation should hold for both vacuum and 
charged configurations. 



3 Extremal solutions: the near-horizon geometry 
3.1 The problem 

The discussions in the previous Section concern the case of nonextremal static black holes possessing a 
nonzero Hawking temperature. However, these findings cannot be used to argue against the existence of 
static extremal configurations. In contrast, static extremal configurations would have zero temperature and 
therefore would not Hawking radiate. Thus the arguments put forward by Emparan et al. |11| and Tanaka 
|12) against the existence of static black holes on the brane would be circumvented. From this point of view, 
arbitrarily large extremal black holes might indeed exist on the brane. 

In principle, the metric ansatz (j2.7l) could still be used in an attempt to numerically construct extremal 
black holes, after replacing the background function F by F(r) = (1 — (r^jr^'^^'^Y ■ (Then the bulk metric 
would still depend on the coordinates r, However, the numerical construction of extremal black holes is a 
highly non-trivial tasll^. Therefore we do not attempt here to construct such bulk configurations numerically. 
Instead, in the following we only address a simplified problem and concentrate on the near-horizon region of 
"potential" extremal black holes. 

Such a restricted study may well provide hints on the properties of the full black hole solutions without 
solving for the full metrico- However, we would like to emphasize that this study could only allow us to rule 
out possible full black hole solutions, while it cannot prove their existence, since for a given near-horizon 
geometry there need not be a corresponding full black hole solution. 

Let us start by recalling that in GR an extremal RN solution in D-dimensions (where later D = d — 1) has 
a near-horizon geometry ^^5*2 x S^~^, which is also a solution of the equations of motion. For a Lagrangian 
density L — R — F^, the corresponding near- horizon line element reads 

ds^ — L\^Qj^^dYj\ + -t'2(Gfl)^^i)-2i (3-1) 

numerical solutions are known to exist, which do not possess a static limit 1251 1261 . 

^^We note that the observed oscillatory behaviour of the configurations is not predicted by perturbation theory (see [7] and 
references therein). However, we suspect that this is a deficiency of the perturbation theory. A somewhat similar situation is 
encountered in Einstein gravity coupled to nonabelian matter fields. There perturbation theory predicted the existence of d = 4 
rotating Einstein- Yang-Mills solitons |27l . Such configurations were, however, ruled out by nonperturbative arguments. 

^^In fact, we are not aware of any successful numerical construction of extremal black holes, obtained as solutions of second 
order partial differential equations. 

^^Such an approach has been used recently to investigate properties of higher dimensional extremal black holes with a 
nonsphcrical horizon topology, see e.g. |28l . 
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where 

-1 _^\J^-'^)^2 t2 _^\iy-0) 2 



T'2 _2(D-3) 2 r2 _ 2(1? -3)3 



dYj2 ~ dx^ /x^ — x^dt^ is the Hne element of the two-dimensional AdS space, and the gauge field is F^-t — Q 



(Note, that the electric charge Q of the bulk RN black holes is Q = QLEq^^/L^/qj^) (up to a volume 



factor).) The entropy of this solution is 

Sd{gr) = '^§^^2{GRy (3-3) 



Turning now to black holes on the brane, we expect that the near-horizon geometry of static extremal 
black holes on the brane in the RS model will - analogously - have a near-horizon geometry ^^5*2 x S"*""^ 
(where AdS2 and S'^~^ have constant radii), which should then be the induced metric on the brane. Thus 
we assume that this metric has a form like eq. (|3.1I) . with, in general, different values for the constant 
coefficients, but recovering eqs. p.2p and p.3p in a certain limit. 

These considerations correspond in fact to the approach of Kaus and Reall [19] used to study d = 5 
extremal black holes, charged with respect to a Maxwell field on the brane. Kaus and Reall [19] found that 
the GR results are recovered for large black holes. (We further note, that Suzuki et al. examined the 
case of extremal black holes in a braneworld with a cosmological constant.) 

3.2 The bulk: equations and asymptotics 

The considerations above lead us to propose the following line element for the near-horizon limit of an 
extremal bulk black hole, 

rfs' - + aiOdY^l + e'rf^iLs, (3.4) 

containing two metric functions a(^) and f{S). Here the coordinate ^ > is the coordinate normal to the 
brane, and the brane is located at some > 0. Thus the coordinate ^ is proportional to sinx in the bulk 
parametrization eq. (|2.7|f^ . 

The functions a(^) and /(^) are solutions of the differential equations 

„ a'^ {d~5){d-4) 1 (d-4)af jd - 4)a' a' f 1 Aa 

4Af 2(d-4), , 2 fa' 
d — 2 4 a 



together with the constraint equation 



,2 2(d-3)(<i-4)a2 1 4a(l-hAa) 4(d-3)aa' 



(3.6) 



These equations have the following exact solutions. 



t2 



/(e) = i + |2' aiO = e + i^ (3.7) 

which corresponds to AdSd in coordinates adapted to a foliation by AdS2 x S"*^^ hypersurfaces, and 

/(O^l + I^U, <.«)^^, (3.8) 



^®We employed a similar metric ansatz and the same numerical methods in our previous studies 1301 131| of asymptotically 
AdS solitons with a nonstandard asymptotic structure. 

'^''The coordinate r in eq. 112.71 1 becomes the coordinate x in the AdS2 parametrization = dx^/x^ — x^dt^ after taking 
the near-horizon limit. 
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which corresponds to AdS2 x H'^^^. 

For the general solutions of these equations smoothness at ^ = requires that a(^) and /(^) possess a 
Taylor series expansion there, consisting of even powers of ^ only, with a(0) > and /(O) = 1. To order 
the small ^ expansion reads 



d-2 + 2aoA ^, ^ (d-2 + 2aoA)((d-l)(d~2) + 2aoA) ^4 ^ 



^ 2(rf-2-ao(d-4)A) 2 2(rf - 2 + 2aoA)((rf - - 2) + 2aoA) ^4 
•'^^^ ao(d-2)2(d-3) ^ a2(^_2)4(d-3)d ^ ' 

One can also construct an approximate solution for large values of the ratio (,/£ (although we shall see 
that configurations with this asymptotics exist for a restricted set of initial data, a(0) > /{d — 1) only), 
which is useful in our analytical study of large black holes. For even d, this solution reads 



aio = ui'\^^+ j2 Mjr-Mqr-'j+o{i/^^'-'), (3.10) 



where the parameter U corresponds to the asymptotic ratio of the radius of AdS2 to that of Also, 
and fk are constants depending on the spacetime dimension d and U only. Specifically, one finds 

_ l + {d-A)U (d-4)(l-t/)(l + (d-4)[/) 

U{d-3) ' (d-2)(d-3)2(d-5)C/2 ' ^^-''^ 

(d - 4)(1 - C/)(l + (d - 4)C/) (2((i - 5)d - 4 + (d - 4)(26 + (3d - 2S)d)U) 

02 — 



3(d - 2)2(d - 3)3(d - 5)(d - 7)t/3 



and 



^ 2+(d-l)(d-4)[/ ^ 2(d-4)(l-t/)(l + (d-4)t/) 
^° C/(d-2)(d-3) ' (d-2)(d-3)2(d-5)t/2 ' > 

_ 2(d - 4)(1 - U){1 + {d- A)U) {d{d - 7) + 8 + (d - 4)(11 + (d - 8)d)[/) 
" (d-2)2(d-3)3(d-5)(d-7)t/3 ' 

their expressions becoming more complicated for higher order fc, without exhibiting a general pattern. The 
corresponding expansion for odd values of the spacetime dimension is more complicated, with log(^/£) terms 
in the asymptotic expression 

/ (rf^/2 £ \ losf 

a(0 = C/|^+ E «'^(^)''+9lo4^(p'-' + (-A^ + /5)(^rM +0(p4), (3.13) 

/(O = 1+ E /.(^)^'= + 2glog(|)(|)'^-^-2M(^)^-3 + 0(^), 

where Ofe and fk are still given by eqs. (13. lip and p.l2p . while q and /3 are two new constants depending on 
£ and d, that can be expressed in a compact form as 

^ ^''^ 3200C/3 ' 42674688 9,d ^ , \ i 

_l-L/2 3(1 - /7)(1 + 3;7)(2 + 3;7) 5(1 - C/)(l + 5J7)(171 + llOOC/ + 1375C/2) 

" T2C7^ ''"^ 800CP + 1778112C/4 '^^^'^ + ' ' ' ' 
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The parameter M in eqs. p.lOp and p.l3|) is a constant, that can be fixed by numerical calcuations, but 
whose value is not of interest in the present context. 

We remark, however, that these bulk solutions are interesting in themselves, since they provide gravita- 
tional duals for some CFTs in a fixed AdS2 x S'^^'^ background given by 

ds^ = e^{u^dT,l + dnl_^). (3.15) 

These configurations have a well defined mass and action which can be computed by using the boundary 
counterterm prescription of Balasubramanian and Kraus |32| . (In fact, both the mass and action are essen- 
tially fixed by the parameter M in the asymptotic expansion, eqs. (I3.10p and p. 131) .) The stress tensor of 
the dual CFT can also be calculated by the method of Myers [33] . A detailed study of these aspects will be 
presented elsewhere in a more general context. 



3.3 The brane 

We now assume that the brane is located at ^ = ,^o, and, following the RS construction, we keep the region 
< ^ < 1^0 of the bulk. The induced metric on the brane is a product of AdS2 and S"'"^, with 

da^ = Ljd^l + Lldnl_^, (3.16) 

where 



Li = vMaO, L2 = Co. (3.17) 

Similar to the case of Einstein-Maxwell gravity discussed above, we shall take a purely electric field with 
Fxt — Q- (Note that, similar to the GR case, the electric charge of the full solution on the brane would be 
only proportional to Q.) 

The Israel junction conditions on the brane, eq. (12. 5p . yield the relation 



a'(Co) _2 (2{d-2) 2d - 7' 



(3.18) 



together with 

= \{d-2){d-2f^ {ym- 1) . (3.19) 

Then the constraint equation p.6p implies that 

1 (d-3)(d-4) _ Q2 , 7d-23 eQ''\ ^g) 



Ll 2Ll L\\d-?, 2(d-2)(d-3)2 L\ 

(which has been used to test the accuracy of the numerical results). Moreover, this relation implies L2 > Li, 
thus the S''*^'^ radius is always greater than the AdS2 radius. 

The entropy of the bulk solution is taken to be one quarter of the event horizon area (note the factor of 
two which originates from the Z2 symmetry of the problem) 

T/, , rio fd-3 



3.4 The solutions 
3.4.1 General results 

We have solved the bulk equations (|3.5p by imposing the initial conditions (|3.9p together with the boundary 
conditions on the brane p.l8|) . p.l9|) for all dimensions between five and nine. 
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Figure 3. The functions a{£) and f{£) are shown for d = 6 and i = 1. The initial parameter assumes the values 
a(0) = 0.4, 0.3, 0.24, 0.2, 0.15, 0.1, 0.05 and 0.02 from top to bottom (left figure) and from left to right (right figure). 

Thus such solutions are likely to exist for any higher dimension d. Different from the strategy employed in 
Section 2, when trying to solve for the full configurations, we have here fixed the AdS length scale £ = 1 
and varied instead the position of the brane, i.e., the size of the black hole as given by the parameter L2, 
eq. (I3T7| . 

The near-horizon region of the extremal black holes has then been constructed in several steps as follows. 
First, we have solved numerically the Einstein equations (j3.5l) by employing a standard ordinary differential 
equation solver. In particular, for a given value ag, we have evaluated the initial conditions (|3.9p at ^ = 10^^ 
and allowed for a global tolerance of 10"^^, when integrating towards large values of (We note, that we 
have not encountered any problems in the numerical integration of the configurations in this Section.) Given 
such a solution of the bulk Einstein equations, in the second step we have used the junction condition (|3.18l) 
to evaluate the position ^0 of the brane. In the final step, the corresponding value of the charge parameter 
Q has been obtained from condition (|3.19p . 

Our results indicate the existence of a single parameter family of solutions of the full problem (i.e., bulk 
plus brane), conveniently labelled by the value a(0). The basic features of these solutions are independent 
of the dimension d. First, for < a(0) < i'^/{d — 1), the bulk solutions do not approach the asymptotic 
form p.lOp . p.l3p . Instead, the function a(^) decreases monotonically and vanishes at some finite radius 
^max, which is a curvature singularity, as seen e.g. by evaluating the Kretschmann scalar. However, such 
configurations are also relevant in the present context, since the junction conditions p.lSp possess a solution 
with ^0 < Cmax- (Thus the singularity is outside the physical manifold.) 

Solutions of the Einstein equations which for sufficiently large values of ^ approach the asymptotic form 
p.lOp and (I3.13|) are present for a(0) > P/{d— 1). In this case, there occurs a maximal value of a(0) as well, 
since the junction condition p.lSp can be satisfied only for a(0) < where the critical value Oc depends on 
the dimension. While Qc ~ 1 ioT d — 5, approximate values for Uc are 0.43, 0.28 and 0.21 for d = 6, 7 and 8, 
respectively. Also, the position of the brane as given by ^0 is a monotonic function of a(0), with both ^o/£ 
and Q diverging as a(0) — Uc, while the parameter U stays finite in this limit. 

The functions a(^) and /(^) of d = 6 bulk solutions are exhibited in Figure 3 for several values of the 
initial parameter a(0). Note that the picture here is generic (as seen e.g. by comparing with the figures 
obtained by Kaus and Reall [19] for d = 5). The exact solution p.8p clearly separates the two different types 
of configurations. 

In Figure 4 we exhibit the parameters Li, L2 and Sd for d = 6 and 7 solutions, normalized with respect to 
the corresponding Einstein-Maxwell results, eqs. p.2p and p.3p with D = d — 1. For any d, the parameters 
ii, L2 and Sd vanish as Q — > 0. 



15 




Figure 4. The ratios Sd/Sd-i, ii/^i{Gfl) a-nd L2/ L2(gr) (shifted by —1) are shown for d = 6 and 7 with 1 = 1. 



However, both Li/Lk^qji) and L2/L2[gr) diverge in this hmit; at the same time, the ratio Sd/Sd-i ap- 
proaches a constant nonzero value {e.g. Sd/Sd-i — 0.14, 0.078 and 0.095 for d = 5, 6 and 7, respectively). 
Moreover, for sufficiently large values of Q, the GR results are recovered. 

For d = 5, we have the choice to consider an electric or a magnetic charge on the brane. However, similar 
to the case of Einstein-Maxwell gravity, one can show that the properties of the solution are the same in 
both cases and the electric-magnetic U(l) duality still hold^. Our numerical results for d = 5 are in good 
agreement with those in [19j . (We note, however, that in [19^ a different parametrization of the metric was 
employed.) 



3.4.2 The large black holes limit 

The observation that ^0/^ diverges as a(0) — Oc, implies that one can use the asymptotic expressions (jS.lOl) 
and (|3.13p to derive a number of simple analytic results for large black holes, and thus to understand the 
large-Q behaviour, demonstrated in Figure 5. Following Kaus and Reall [19 , one starts by assuming an 
asymptotic expression for the parameter U (which enters the large-^ expansion (|3.10p . p.l3|) ) in terms of a 
power series in e = i/£,o- Then the first junction condition (I3.18P implies that 

U 1 (d-3)(4d-ll) 

^ " (d-4)2 4(d-2)(d-4)2' ■ ^-^-^^^ 

This shows that the metric function a(^o) (*-e., the size of the AdS2 part of the brane metric) has the 
following approximate form on the brane 

« - ^1 - Jd^^ - 4(d-2)(d-4)2 + • ■ • • ^^■'^^> 
The second junction condition p.lSp gives for the charge parameter the expression 



d-3 t {d-3){7d~18) 



2(d-4)3e2 8(d-2)(d-4)3 



r + . (3.24) 



^®We thank H. Reall for a correction on this issue. 
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Inversion of eqs. (|3.22p and p.23p yields e as a function of . One obtains 



, 2(d-4) 6d-19 , 

+4(d-2)(d-4)2^ ^^-^^^ 
,_ 2(d-4)^ 7d-18 , 

^2 - (d-3) ^ + 4(d^^ + ■ • ■ ' ^^-^^^ 

where comparison with p.ip (with D — d~X) shows agreement with GR in the hmit of large Qjt. (Note 
that the leading order corrections are strictly positive.) 

The same holds also for the entropy of these solutions, which, in d > 5 dimensions, is given by 

a _g^rf-3 t ,d-3f. 3(d - 3) (d - 4) \ 

^'-^^Jd^f'^ V 2(d-2)(d-5)ijj+--- • ^^■^^> 

The first term here corresponds to the entropy of an extremal RN black hole in c? — 1 dimensions. 
(Recall that Gd — -^^Gd-i-) With help of relation (I3.24p . one can reexpress Sd as a function of the charge 
parameter Q. As seen from the large-Q region of Figure 4, the analytical and numerical results are in 
excellent agreement. 

For d = 5, the leading order correction for the entropy as a function of Q contains a log term, which 
is absent in other dimensions. The corresponding relations are given in |19j . together with some relevant 
plottB. 



4 Further remarks 

With this work we have addressed the issue of black hole solutions in the RS infinite braneworld scenario, 
allowing the black holes to carry electric charge on the brane. In the first part of the paper we have considered 
vacuum black holes and nonextremal charged black holes on the brane. Employing a different numerical 
technique by solving the set of elliptic partial differential equations in the full bulk, ranging from the brane 
to the AdS horizon, we have obtained results that fully support the claim by Yoshino |^ for the nonexistence 
of static vacuum black holes on an asymptotically flat brane in the RS infinite braneworld model. This 
conclusion does not change when (nonextremal) solutions are considered, which are charged with respect to 
a Maxwell field living on the brane. Although "approximate" solutions appear to exist for sufficiently small 
brane tension, these configurations are very likely only numerical artifacts. 

One should emphasize that, of course, the numerical results cannot be used to prove the nonexistence 
of static black hole solutions on the brane. To clarify the issue of the existence of black holes on the 
brane one would thus need to either find a full analytic solution or, alternatively, to provide a rigorous 
theoretical argument for the absence of such black hole solutions. Lacking both, however, we think that a 
natural interpretation of our and previous numerical results is provided by the conjecture that nonextremal 
braneworld black holes would necessarily be time-dependent [TTl [T^] . 

The situation is different for extremal black holes, since the conjecture put forward by Emparan et al. [TT] 
and Tanaka [H] does not forbid the existence of static extremal black holes. Thus, in principle, the presence 
of a second global charge, apart from the mass, could allow for the existence of extremal black holes also 
in a braneworld context. While the numerical construction of localized extremal black holes still represents 
a numerical challenge to be met, we here have considered the less ambitious task of constructing only the 
near-horizon geometry of extremal solutions with a Maxwell field on the brane. In particular, we have found 
that the GR predictions in d — 1 dimensions are reproduced for extremal black holes which are sufficiently 
large as compared to the AdS scale. 

One should mention, though, that finding local solutions in the vicinity of the horizon does not guarantee 
the existence of the corresponding global solutions. (Chen et al. [34 , for instance, give an example where the 
physically relevant global solutions are absent despite the presence of a closed form near- horizon solution.) 

^^For d = 5, the electric charge is Q = -pfQ- 

^1 
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In our opinion, any progress in this direction would require the development of a consistent numerical scheme 
capable to achieve the explicit construction of the bulk extremal black holes. 

However, another physically interesting situation to consider in the context of the RS infinite braneworld 
model concerns the case of static localized solitons on the brane. With no event horizon present, solitons do 
not possess intrinsic thermodynamical properties. Thus the conjecture of Emparan et al. and Tanaka 
[T^ has no obvious bearing on the existence of solitons on the brane. 

Let us thus consider as the simplest example of (possible) solitons on the brane in the RS model, the 
case of brane world Q-balls based on the matter field Lagrange density 

-Cq = a^$*9^$ + C/(|$|). (4.1) 

The matter field in this case is a complex scalar field $ with a harmonic time dependence and a non- 
renormalizable self-interaction described by the potential J7(|<i>|). The fiat spacetime solutions of this theory 
were considered for the first time by Coleman 35 , while their Einstein gravity generalizations were discussed 
in [36l [37l |38]. These gravitating nontopological solitons describe localized particlelike objects with finite 
energy. 

We have attempted to construct d — b spherically symmetric gravitating Q-balls in the RS model, where 
the scalar field is confined to the brane and described by the simple ansatz 

$ = f{r)e-''"\ (4.2) 

where we have followed a basically similar approach for these gravitating Q-balls as for the static black hole 
solutions. In particular, we have tried to solve the bulk Einstein equations numerically for the metric ansatz 
(|2.7|) with ro = 0, where the boundary conditions at r = 0, r = oo and x = are still given by p.20p - ()2.22|) 
(with 7'o = 0) and (I2.5p . respectively, and where the energy-momentum tensor is given by 

t^. = + - g^.Cg. (4.3) 

However, unexpectedly, our attempts to construct gravitating Q-balls on the brane have not been suc- 
cessful. As in our attempts to construct static black hole solutions on the brane, we have failed here to obtain 
reliable numerical gravitating Q-ball solutions on the brane. Interestingly, the reason for the numerical in- 
accuracy of the configurations and the lack of the numerical convergence again resides in the unnatural far 
field behaviour of the metric functions close to the AdS horizon. In fact, this unnatural far field behaviour 
is completely analogous to one observed for the configurations supposed to describe static black holes. 

Clearly, this behaviour is hard to understand, since there is no clear a priori reason for solitons not to 
exist on the brane. However, if we interpret our results for the static (nonextremal) black hole solutions as 
implying, that there are no such solutions on the brane (satisfying the given set of equations and boundary 
conditions), we must - consistently - draw the same conclusion for the gravitating Q-ball solutions. While it 
might be interesting to attempt the construction of other types of solitons on the brane, the present results 
strongly discourage such a construction, since the same problems in the vicinity of the AdS horizon are likely 
to arise. 

How do we then judge the issue of black holes and regular solutions on the brane? There are still the 
positive results from the near-horizon black hole solutions in the extremal case, but these extremal near- 
horizon solutions avoid the problematic region in the vicinity of the AdS horizon. Thus a full calculation of 
extremal solutions is clearly called for, where one attempts to smoothly extend the near-horizon geometry 
of the extremal solutions into the asymptotic region, in spite of the tremendous numerical challenge. The 
final outcome might, however, be that one is led to conclude, that there are no extremal black hole solutions 
on the brane, either. 
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A The numerical method 



The set of three coupled non-hnear ehiptic partial differential equations for the functions A, B and C has 
been solved numerically for the (p, x) coordinate system introduced in Section 2.1.1, subject to the boundary 
conditions t^nU\ - (^n^ . 

The first step has been to introduce a new radial variable x = p/{l + p) which maps the semi infinite 
region [0,cxi) to the finite region [0, 1]. This involves the following substitutions in the differential equations 

pu^p — >{l-x)f^^, p^u^pp — > {l-x)'^u^x^-2{l-x)u^^ (A.l) 

where u denotes any of the unknown functions A, _B, C . 

Next the equations for these functions are discretized on a non-equidistant grid in x and x- Here we have 
considered various grid choices, the number of grid points ranging between 180 x 30 and 80 x 70. The grid 
covers the integration region < x < 1 and < X < 7r/2. Typically, the mesh in the a;— direction has been 
denser in the near-horizon region [x = 0) and for values of x close to the AdS horizon. Most of the x-meshes 
employed have been equidistant. 

All numerical calculations have been performed by using the programs FIDISOL/CADSOL J23]. Here we 
shall briefly review its basic aspects. The code requests the system of nonlinear partial differential equations 
in the form 

P{x, y, U, Ux, Uy, Uxy, Uxx7 Uyy) — 0, [A..2) 

subject to a set of boundary conditions on a rectangular domain. (For convenience, we have used the 
notation x = U-) Besides the set of equations, FIDISOL/CADSOL requires the boundary conditions, the 
Jacobian matrices for the equations and the boundary conditions, as well as some initial guess for the 
functions. The Jacobian matrices are generated by simple differentiation of each equation P with respect to 
and Uyy. 

FIDISOL/CADSOL uses a Newton- Raphson method. The numerical procedure works as follows: for an 
approximate solution u^^\ P{u^^^) does not vanish. The next step is then to consider an improved solution 

= +iuAm, (A.3) 

supposmg that P(w(i) + wAu) = (with w a relaxation factor, which is usually chosen as w — 1). The 
expansion in the small parameter Au gives to first order 

f)P 

= P(u(i' + Au) « F(m(i)) + ^(m(i))Au + . . . . (A.4) 

ou 

This equation is then used to determine the correction Au^-^^ = Au. Repeating these calculations iteratively 
(y(3) — y(2) _|_ Ay etc.), the approximate solutions will converge, provided the initial guess is close enough to 
the true solution. In each step, a linear system of equations is solved, and the residual ||P(m^*^)|| decreases 
by a factor of approximately 10 — 20. The iteration stops after i steps, when the Newton residual P(u'^'^) is 
smaller than a prescribed tolerance. Clearly, it is essential to have a good first guess, to start the iteration 
procedure. 

The package FIDISOL/CADSOL provides also error estimates for each function, which allows to judge 
the quality of the computed solution. The errors are computed on the "consistency level" , namely, the 
discretized Newton residual, and as discretization error terms in x, y directions. The discretization error is 
estimated through the difference of difference quotients. For example, in (jA.4[) . the derivative of the solution 
u and of the correction function Au are discretized by a difference method with arbitrary consistency 
orders. Derivatives are replaced, for example, in the form Uxx Uxx.d + dxxi Auxx Auxx,d, where 
the index d means "discretized". dxx is the estimate for the discretization (or truncation) error of Uxx, 
defined as dxx — Uxx,d.next — Uxx,d, where the index "next" denotes the next higher member of the family of 
nonequidistant backward difference formulas. The discretized Newton residual decreases with the number 
of Newton- Raphson iterations, while the discretization error terms depend on the grid size and the used 
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consistency order, i.e. on the order of the discretisation of derivatives (in our work, this order was six). 
Furthermore, the error terms are used for the determination of stopping criteria for the Newton-Raphson 
method. Further details on the numerical method and explicit examples are provided in |23) . 

In this scheme, the input parameters are the event horizon radius tq and the value i of the AdS length 
scale which form the dimensionless parameter L = r^ji. In our approach we set rp = 1 and we start with 
the Schwarzschild-Tangherlini solution as initial guess (i.e., L = and A — B — C = Q^. Then we increase 
the value of L slowly. The iterations converge, and, in principle, repeating the procedure we obtain in this 
way solutions for higher values of L. In some of the calculations, we interpolate the resulting configurations 
on points between the chosen grid points, and then use these for a new guess on a finer grid. 
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